Finite-temperature correlations in the trapped Bose gas 
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j_j ' Previous functional integral methods for translationally invariant systems have been extended to 

O . the case of a confining trap potential. Essentially all finite-temperature properties of the repulsive 

^j ■ Bose gas in a paraboloidal trap can be determined this way. New analytical results reported here are 

for the finite-temperature two-point correlation functions below the critical temperatures in d = 3, 2 
and 1 dimensions. Only for d = 3 are correlations both long range and coherent - at long range 
consistent with the existence of a coherent state of the condensate in the trap. The two-dimensional 
^ I condensate is marginally stable in the sense that correlations decay by a power law. 

^ ! 03.75.Fi, 05.30.Jp 

O ' The observation of Bose condensation in vapours of alkali atoms [El-pl held in magnetic traps, and recently Hi 

a , in atomic hydrogen, has stimulated enormous interest, both experimentally and theoretically, in this phenomenon. 
The coherence properties of condensates in such traps are of exceptional interest especially perhaps because of the 
4_> I opportunity now presented for realisation of atom lasers [^-|7| . Much of the theoretical work has been concerned with 
^ ' solutions of the c-number Gross-Pitaevskii (GP) equation in the presence of the paraboloidal potential describing the 
c/3 \ trap p-[lO|. Without the trap this equation is also called the Nonlinear Schrodinger (NLS) equation |l^ which, in 
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one space dimension {d — 1), can be solved exactly at both c-number and quantum level, the latter including [02 13 
calculation of the finite-temperature correlation functions for the repulsive case (coupling constant g > 0). For d = 1 
as well as for the higher dimensions, the finite-temperature properties of the quantum NLS equation have otherwise 
been extensively analysed [|^ by functional integral methods. This way it is established e.g. that without any trap 
2 I no long-range correlations arise for d < 3. But, on the other hand, for d = 3 correlations are both long-range and 
Q ■ consistent with the condensate being in a coherent state |13] . 

O I Thus far it has remained an open question whether the presence of a trap potential will induce long-range order for 

d < 3 also, while, despite P,0|, coherent-state behaviour for d = 3 in traps had not been demonstrated theoretically 

I until the very recent work of [ |l6[ . In [|l6| first and second order coherence functions are calculated from the GP 

^ , equation in the form of series expansion. In this Letter we give an ab initio approximate calculation of the first order 

^^ ' coherence function for d = 3,2, and 1, giving explicit analytical formulae in each case. We show that, even in the 

C^ i presence of a paraboloidal trap, there is long-range order only for d — 3; and we indicate how, and only for d = 3, 

\l ' the condensate in the trap may, but only asymptotically, be approaching a coherent state. To these ends we have 

^•— ^ [ extended the previous |14| , |l7| functional integral methods to the case when a trap is included. An important new 

^^ ■ aspect is then that the trap breaks translational invariance; and this introduces wholly new many-body theoretical 

Qv ' problems. Thus, rather than using periodic boundary conditions (b.c.s), we must expect to impose vanishing b.c.s at 

"•-••, infinity, and so to work at 'zero density' |l8[ . In this Letter we simply divide ||l9,pQ| the total field ip of the system into 

C^ two parts: ip = ipo -^'>Pi- Here ipo becomes for T < Tc the condensate field and -01 is that due to thermal fluctuations. 

C Both fields are assumed to vanish at infinity in order to have them largely confined in the region of the trap. The 

I functional integral method still provides a framework within which, in principle, all thermodynamic properties of the 

■ O _ trapped bose gas can be determined. So far we have obtained, e.g., expressions for Tc for each oi d = 3,2, 1. At 

5 ■ mean-field level results at lowest order in the coupling constant agree with those reported ||] already, but they can be 

extended to higher order. However, the main result reported in this Letter is the calculation of the finite-temperature 

two-point correlation functions for each of d = 3, 2, 1. Because there is no translational invariance these correlation 

functions no longer depend solely on the difference of two position vectors. It can still be concluded that long-range 

order arises, for T < T^ in d = 3, and the first-order coherence function asymptotically approaches unity. In d = 2 

5J] ■ and for T < T^ the condensate is marginally stable for correlations decay algebraicly, namely as a power law. In d = 1 

Ci I correlations decay exponentially for T < Tc, and we have not yet analysed any T = limit. For T > Tc there is a 

Gaussian decay in all dimensions. 

Our calculational procedure is such that we first draw a correspondence between the quantum fields ip{r, t), V'^(r, t) 
for Bose particles and the c-number fields ip{r,T),ip{r,T) (r is the Wick rotated time t -^ — ir, and varies be- 
tween and P = (kBT)~^ ). In this way we can write the finite-temperature correlation function G(ri,r2) = 
(TT--0(ri, ri)V'^(r2, T2)) (where Tt- means a thermal ordering in r) as the ratio of two functional integrals, 

G(ri, r2) = Z-^ [ e^V(ri, ri)V;(r2, T2)i?V^V5, (1) 



in which Z is a. partition function Z = f e^DipDil;. The chosen c-number action S is 

S^j dT d'^ri.'iJj{r,T)KiP{r,T)-^i>{r,T)'iP{r,T)^{r,T)iP{r,T)\ . (2) 

The boundary conditions are vanishing at infinity for r and periodic, period /3, for r. The action S yields the 
quantum many-body problem at T > for a gas with repulsive pairwise (5- function interactions of strength g in R"*. 
The differential operator K — dr ~ H, and H — — ^^V^ + V{r) — /x; /i is the chemical potential and V{r) = ^O^r^ 
is the paraboloidal trap potential, taken with spherical symmetry for simplicity. 

As explained we also simplify by setting ip{r, r) = ipoii^) + V'i(r, t) for c-number fields ijj, ipo, ipi and likewise for tp : 
ipoi^) will not depend on r. In this Letter we shall only consider terms in S upto quadratic (bilinear) in V'liV'i ^nd 
can therefore explicitly integrate out the thermal fluctuations. This way we arrive at Seff[ipo, ipo] = In/ e^DiJjiD^pi, 
and 

Seff [ipo, V'o] + PPnaip) 

d''r{4>o{r)[—V^ + ^l- V{v)]i,,{v) - Uo{v)Mv)M^)i,o(v)}. (3) 
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In Seff, V{v) = V{y) + 2gpnc{'r) is a renormalised trap potential, while p„c(r) is the density profile of the thermal 
particles in the ideal gas approximation - as is consistent with terms only quadratic in ipijipi retained. More precisely, 
at this level of approximation (Hartree-Fock-Bogoliubov, HFB) this density profile derives from the fundamental 
solution of the d+ 1-dimensional operator K: KQ{r,T; r' ,t') — —S{r— r')6{T ~ r'), and the thermal Green's function 
Q can be expressed in the form 

e(r, r; r', r') = ^ !^^g^e^"(-^'), r > r'. (4) 

n 

For d — 3 the vectors n — (ni,n2,n3), and the Mn(r) and E„ are the eigenfunctions and eigenenergies, respectively, 
of the d = 3 harmonic oscillator Hamiltonian H; and similarly for d = 2,d = I. Then /0„c(r) = Q'{r,T;r,T), 
where prime means n = (0,0,0) is omitted. The free energy of the thermal particles Fnc{p) is simply FncifJ-) — 

At this point we are already able to calculate the critical temperatures Tc- The leading i.e. zeroth order term is 
found by replacing /9„c(r) by a constant p„c(0): this defines the renormalised chemical potential A = /i — 2gpnci0), 
and A = determines T^. We find this way all of the zeroth order expressions in d = 3, 2, 1 as given e.g. in ||^. At 
first order in g we can add the appropriate terms arising in Seff . Beyond this we also need to include fluctuations ^i 
to an order higher than quadratic. 

The free energy F of the trapped Bose gas is calculated from —(3F{ijl) = In J e^''ffDtpQD'ipQ. By steepest descents 
for large (3 (low T) we find that 

F(/.)=F„,(/i)-|y'd'^r|$(r)|4. (5) 

In Eq.(|^) the fields <&,$ are the quasi-classical fields satisfying the extremum condition 5{Seff[^,^]) = 0. This 
condition is equivalent to the stationary GP equations 

|iv2<l>(r) + (a. - F(r))<i>(r) - g \ $(r) p <i>(r) = 0, (6) 

and the similar equation for $. At this quasi-classical approximation we already find through the presence of ^(r) 
the HFB corrections to the GP equation introduced earher |^. 

We turn next to the calculation of the finite-temperature correlation function Eq. (|l|). By integrating out the 
thermal fluctuations included upto terms quadratic in ipi ,"01 we flnd that 

G(r,.r.). ^'°-7f;]*°;'''f°°'^° .a(r..r,)- (7) 

J e^offDipoDipo 

At low enough temperatures these remaining functional integrals can be evaluated by steepest descents where again 
we work consistently at the HFB level. At this level the correlation functions can be expressed in the form 



C(ri,r2) ~ g-S^ff['i'o,'S'o] + S^ff[<S>i,<S>i]+\n<i>i(ri)'S>i(r2)^ rg\ 

and the fields $o, ^o satisfy the two stationary GP equations exemphfied by Eq. (g). Evidently the fields $1, $1 are 
determined by S{Seff[^i, $1] + ln$i(ri)i'i(r2)) = 0, and this variational equation also leads to a pair of GP type of 
equations with however additional sources from the ln($i$i). This pair of equations is 



-lv^^,(r>- 


~{^i~V{r))^l{r)+g<^l{r)^l{r) 


-lv'^,(r>- 


-{^i-Vir))^l{r)+gmr)^l{r) 



(9) 
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/3$i(r2) ' 
5{r-ri) 

/?$i(ri) ■ 

In this Letter, for simplicity, we shall solve both these equations Eqs. (^ and the Eq. (^) only at the Thomas-Fermi 
(TF) approximation: this is expected to be valid at low enough temperatures |8|-[l0[. The stationary GP equation 
Eq. dg) is thus reduced to a simple algebraic equation and from it we then easily find the expected inverted paraboloidal 
density profile which is P,pi[ 

po(r) = I'oW^oW = -(/i - V{r))e{^I - V{r)), (10) 

9 

in which 6 is the Heaviside step function. Evidently /9o(r) can be interpreted as the condensate density expressed 
in terms of order parameters $0,^0 ■ The radius of the condensate Re can now be determined from the condition 

The solution of Eqs. (^ is more complicated. Notice first that the fields $1, $1, appearing in these equations will 
also be complex valued fields in general. But they are two independent fields with independent variations. We can 
therefore seek first of all solutions in the form <I>i(r) = ^p(r)e'^'''\ $i(r) = ^p(r)e~'^^'"\ where (/)(r) can be complex 
valued, but will be found below to describe the real contribution to the correlation functions of the complex phases 
of the wave functions. We can furthermore assume that, away from the boundaries, p{r) is a slowly varying function 
of position r so that V^.^ and Vy^ are both small and can be neglected. This will not be true of p(V^)^ or pV'^cj), 
so that Eqs. (ph can be taken as 

5p(r)-(^-f(r))-^(V0(r))2=O (11) 

^V'Hr) - ^^JT^Mr - n) - -T^'5(r - r^). (12) 

2m 2/3p(ri) 2/3/)(r2) 



The first of these equations Eq. ( |ll| ) has the solution p(r) = po{r) + ■^^{'^4') ■ Within the TF approximation, p(r) 
in Eq. ( p^ ) is then po{ic)- We then express the solution of this equation in terms of a function /(r,r'): (/)(r;ri,r2) — 
/(r, ri) — /(r, r2). The functional form of /(r, r') depends on the dimensionality of the system. We find that 

fir,r') = -^R {d^l) (15) 

with a = -^ and i? = | r — r' | . It is already clear that the correlation functions can no longer depend on i? = | r — r' | 
alone: they depend also on both of ri and Y2 separately, consistent with the breakdown of translational invariance 
induced by the trap. 

We consider first the correlation function in d = 3. In this case the points r = ri and r — Y2 in 0(r;ri,r2) are 
singular and introduce a divergence problem |22|] . This difficulty can be avoided first of all by considering a first-order 
'coherence function' G^^^(ri,r2) which we define here (and compare e.g. |lq,E3]) as 



where C was defined at Eq. dzh. For this way G'-^-' is both finite and well defined because we find identically the 
same singularities appear [Q in the direct calculation of the order parameters {^{ri,Ti)), {ip^r2,T2)). Notice that 
in Eq. ( p^ ) we have already replaced {tp), the order parameter of the trapped Bosegas by (ipo) since the average 
over thermal fluctuations vanishes: {tpi) = 0. Notice too that the solution of Eq. (|lj) is only determined within a 
harmonic function, and that the only harmonic function which satisfies the b.c.s. and is non-singular is a (complex) 
number. This constant is identically cancelled in G^^' but introduces an undetermined phase, as expected [p3-E6|, to 
the order parameters. For T < Tc, when the order parameter is nonzero, we find that 

G(^^(ri,r2) ~ e-~iifiri,r2)+f{r2.r,)) ^ ^in^ jz ^ (^ij) 

where fi^^ = pg"^(ri) + pQ^{r2)- Evidently G'^^^(ri,r2) -^ 1 for large R, thus indicating long-range order and long- 
range coherence, and there are thus these features of a coherent state in this limited sense [ P7[ . The coherence length 
is given by ^^° and depends on both ri and r2 separately. Notice that we have assumed ri and r2 are not close to 

the boundaries of the condensate so that always po{r) > in p^^ . 

For d — 2 the singularity in /(r, r') is logarithmic and the divergence is renormalizable. For d = 1 the function / 
is nonsingular. Thus we can directly evaluate the correlation functions, and find that 



G(ri,r2) ~ VPo(ri)po(r2)e-^^3^i"« (d = 2), (18) 



G(ri,r2) ~ Vpoiri)poir2)e-^m'' (d - 1). (19) 

It is obvious that these correlation functions both vanish for large R and that there is no long-range order in d = 2 
or in d = 1 nor will there be any coherence. In the case d = 2 the condensate is marginally stable in that correlations 
decay algebraically, namely by a power law. The exponent of this power-law is proportional to T so that at very low 
temperatures, correlations may thus prevail over almost macroscopic distances. In real magnetic traps for d = 3 we 
still expect the condensate to be stable even for extremely anisotropic trap potentials, as in e.g. the experiments H] 
on atomic hydrogen. The three correlation functions Eqs. (|lj-|l9|) coincide with those obtained under translational 
invariance without the trap to the extent that for fJ —* 0, V{r) — >0, and we can expect 2g/9„c(r) — *-Const= 2gpnc- 

Thus in summary we have demonstrated that the functional integration techniques developed for translationally 
invariant Bose systems can be extended to Bose gases in a confining trap potential, and form also in this case of a trap 
a convenient framework in which to consider the thermal properties of the condensate. We have shown in particular 
that true long-range order only arises in d = 3 and that the condensate then appears asymptotically to be in the 
coherent state essential for atom lasers. We expect to define multi-point correlation functions similarly to Eq. ( [iq ) 
and find these ~ 1 for large separations of the points. 
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